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Using in-field single crystal neutron diffraction we have determined the magnetic structure of 
TbMnOa in the high field P\\a phase. We unambiguously establish that the ferroelectric polar- 
ization arises from a cycloidal Mn spins ordering, with spins rotating in the ab plane. Our results 
demonstrate directly that the flop of the ferroelectric polarization in TbMn03 with applied magnetic 
field is caused from the flop of the Mn cycloidal plane. 

PACS numbers: 61.12.Ld, 61.10.-i, 75.30.Kz, 75.47.Lx, 75.80.+q 



The antisymmetric Dzyaloshinski-Moriya (DM) 
interaction [TJ [2] between two spins, Si, Si+i separated 
by r^i+i, provides for a natural coupling between 
magnetism and ferroelectricity with the spontaneous fer- 
roelectric polarization given by P s ~ x (Si x Sj+i) 
[3J |U [S]. This mechanism generates ferroelectricity 
in a wide variety of magnets such as i?MnC>3 per- 
ovskites with R=Gd, Dy, and Tb,|3[7], spinel chromate 
CoCr 2 04,[8] spin-chain cuprate LiCu202,[9] and hueb- 
nerite MnWOi [10| ITT] . In the iflVln0 3 manganites the 
DM interaction results in cycloidal order of Mn-spins 
giving a spontaneous ferroelectric polarization along the 
c-axis (P\\c) (Fig. []Ja)). The application of magnetic 
field results in the flop of the polarization from the c- to 
the a-axis and highlights a novel control of one ferroic 
property by an other [B]. It has been assumed that this 
change in the direction of the polarization reflects the 
flop of the Mn spin cycloid, implying that the antisym- 
metric DM interaction continues to be responsible for the 
polarization (Fig. [ljb)) [3j [12] . However, in this high field 
P||a-phase, the commensurate magnetic wave vector 
for Mn is also compatible with other magneto-electric 
mechanisms such as exchange striction[T2] HH [TS]. In 
this letter we present a determination of the magnetic 
structure of TbMnC>3 in a high magnetic field in the 
commensurate P\\a phase. We find that the magnetic 
structure of Mn-spins is characterized by an afr-cycloid 
that accurately describes the direction of the observed 
ferroelectric polarization via the DM interaction. Our 
finding validates the model that the polarization flops 
found in the perovskite manganites result from the flop 
of the Mn-spin cycloid. 

The manganite TbMnC>3 crystallizes in the orthorhom- 
bic perovskites structure Pbnm. On cooling, below 
Tjy=41K Mn spins order incommensurately point along 
the magnetic wave vector t ~ 0.275b* E]. On further 
cooling below Tg=28K, a c-axis component of the Mn 
moment orders with a phase shift of ir/2 with respect to 
the 6-component so as to form a cycloidal structure where 



Mn-spins rotate within the be plane and around the a- 
axis as shown in Fig. [l|a)[5]. The axis of spin rotation 
defines the DM interaction, Si x Sj+i, while the distance 
r^j+i is parallel to the modulation vector t. For this 
type of spin order, inversion symmetry is broken yielding 
for i?=Tb and Dy, a ferroelectric polarization along the 
c-axis as indeed is observed (P s ||c ~ a x b)[51|B]. 

It is tempting to assume that the flop in the ferroelec- 
tric polarization arises from a flop in the Mn-spin spiral, 
but so far there is no experimental prove for this. Fur- 
thermore, the field dependence of the Mn-spin spiral is 
difficult to analyse as single-ion anisotropy terms of Mn 
and R need to be taken into consideration as well as the 
R-Mn interaction. The fact that in this P\\a phase Mn- 
spins order commensurately with r = |b*|15|,H6] render 
the cycloidal flop model even more complex. For such a 
commensurate spin structure it has been proposed that 
spin frustration and super-exchange induce lattice distor- 
tions that break inversion symmetry thereby generating 
the observed direction of the polarization at high field 
[13l [HJ [15]. This exchange striction model can be ap- 
plied only to a commensurate order, ruling out its va- 
lidity for i?=Dy|T7]. The observation that infrared elec- 
tromagnon signals are always polarized along the a— axis 
irrespective if P s is parallel to the a- or c-axis further 
adds to the debate of the high field magnetic phases 
in these manganites [13]. A cycloidal Mn magnetic or- 
dering that yields P\\a in zero field has been observed 
for Gdo.7Tbn.3Mn03[]35], however, there is no direct ev- 
idence that the magnetic field flops the Mn-cycloid to 
yield a P\\ a polarization. It is therefore pressing to estab- 
lish an accurate model of the high field magnetic struc- 
ture of these multiferroic manganites. 

TbMn03 single crystals were obtained by re- 
crystalizing a ceramic rod under Ar atmosphere using an 
optical floating-zone furnace. The field and temperature 
dependence of the magnetic propagation wave vector was 
measured on the E4 double axis neutron diffractometer 
at the BENSC facility of the Helmholtz Zentrum Berlin 
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FIG. 1: (Color online) (a,b) Illustrations of the two cycloidal 
Mn magnetic structures proposed to cause ferroelectricity in 
TbMnOs. In both cases the magnetic propagation vector r 
is parallel to the 6— axis. In zero field (panel (a)) cycloidal 
order, Mn spins rotate around the a-axis (Si x Si+i) and ro- 
tate wholly within the bc-plane. The ferroelectric polarization 
via the antisymmetric DM interaction is produced along the 
c-axis. The high-field magnetic structure determined in this 
work is shown in panel (b), it yields a P||o ferroelectric polar- 
ization that arises from an ab cycloid where Mn spins rotate 
around the c— axis. In both panels we also show the magnetic 
ordering of Tb-spins. In the zero field case (a), the mag- 
netic propagation vectors of Tb and Mn-spins are clamped 
and Tb-spins point along the a-axis forming a SDW [5]. The 
canted antiferromagnetic ordering of Tb-spins for H\\b=5T 
determined in this work is shown in panel (b). (c) Here we 
depict an anharmonic ab-spiral where the phase difference be- 
tween spins 1 and 2 is u) — n/2. Is such a case the angle 
between spins 1 and 2 is different from that between 2 and 
1+b yielding an alternating scalar product along the b— axis. 
Amplitudes of Tb and Mn spins are not to scale. 



using a neutron wave length of A=2.45A and a A/2 fil- 
ter placed in the incident beam. Here magnetic field was 
applied horizontally along the b— axis using the HM1 su- 
perconducting cryo-magnet. Due to the limited view of 
the sample in this magnet, integrated intensities of Bragg 
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FIG. 2: (Color online) Single crystal neutron diffraction mea- 
surements from the E4 diffractometer. Here we show scans 
along (0,k,l) in reciprocal space as a function of temperature 
measured in a 5T field applied along the b-axis. The inten- 
sity of these scans is plotted in color coding with correspond- 
ing scales above each panel, (a) Portion of the data showing 
the temperature dependence of the A-mode reflection (0,e,l). 
Here the wave number e varies from 0.282 at Tn to 0.266 at 15 
K before it locks discontinuously to the commensurate value 
of | below 11K. The weak third harmonic reflection was also 
observed in these scans to disappear at this transition. These 
data are shown in an enhanced color scale on the same panel 
and in the same location in the T — Q map. (b) Temperature 
evolution of the Pbnm forbidden reflection (0,0,1). 



reflections could only be measured within certain regions 
of the Okl lattice plane, a situation that prohibits accu- 
rate analysis of the magnetic structure. To overcome this 
problem we conducted measurements using the D23 neu- 
tron single-crystal diffractometer installed at the Institut 
Lauc-Langevin (ILL) with A = 1.28lA, that is equipped 
with a lifting detector that allows measurements of Bragg 
reflections above and below the scattering plane. Mag- 
netic field was applied using the 6T vertical field super- 
conducting magnet with an asymmetric vertical opening 
angle of -5/+ 10°. For these measurements a single crys- 
tal of TbMnC>3 was cut into a parallelepiped with dimen- 
sion 3. 0x2. 9x3. 7mm with each face perpendicular to one 
crystallographic direction. The crystal was oriented with 
the 5-axis parallel to the field. This geometry allowed us 
to measure reflections with k from -1.8 to 0.25. In to- 
tal 138 independent nuclear Bragg reflections consisting 
of 310 individual reflections were collected at H = OT 
and T=8.5K within a range of (0.07 < sin(6>)/A ^ 0.69). 
At 8.5K the magnetic field was applied along the &-axis 
to 5T. In these conditions 140 commensurate (t = 0) 
independent Bragg reflections consisting of 160 individ- 
ual reflections were measured along with 64 independent 
reflections with r = jb*. Analytical absorption correc- 
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tions for each reflection were made by the Xtal suite of 
program, while analysis of the magnetic intensities was 
performed with the Fullprof code. 

In Fig. [2^a) we show measurements of the (0,e,l) re- 
flection as a function of temperature with H ||6=5T mea- 
sured on the E4 diffractometer. The data shows that Mn- 
spin ordering is first observed at T/v =4 IK with the wave 
number e decreasing on cooling. At T C =11K, e changes 
discontinuously to yield a commensurate wave vector of 
t = |b*[T21 HH1 123- Published polarization data shows 
that a P\\c state develops below T$ while the polariza- 
tion flops to P\\a at T c , coinciding with the transition 
to the commensurate magnetic phase[21 . In the same 
data we observe that the third harmonic reflection (0,1- 
3e,l) rapidly changes its position and disappears also at 
T^Fig.ga)). On cooling below T N we find an enhance- 
ment in the intensity of several nuclear reflections and 
the appearance of forbidden reflections such as the (0,0,1) 
shown in Fig. [2jb) . As we discuss below these effects arise 
from the ordering of Tb-spins with magnetic propagation 
vector T Tb = 0. 

To determine the magnetic structure of TbMnC>3 in 
the P\\a state we utilized the D23 diffractometer. Here 
the sample was cooled in zero field to 8.5K and then 
a H\\b=5T was applied so as to enter the commensu- 
rate P\\a phase that was confirmed by measurements 
of the magnetic wave vector. At this field and tem- 
perature we find that the most intense magnetic reflec- 
tions with wave number e = ~, have extinction condi- 
tion h + k =even, l=odd ( A-mode) |24j . while reflections 
with h + k =odd, I =odd (G-mode) were considerably 
weaker. For space group Pbnm and wave vector r ^5 jb* 
there are four irreducible representations (irreps) T of 
the magnetic symmetry for the Mn-ion [551 The 
A mode reflections are contained only in irreps T\, T2 
and T 3 . Analysis of the data using only A-mode reflec- 
tions and a single irrep did not result in a satisfactory 
fit to the data. This lead us to consider combinations 
of representations. Of the possible combinations only a 
model using Fi ® r 3 produced a satisfactory result with 
i?(F 2 )=7.9% and R W (F 2 )=8.5%. This coupled irrep has 
the form of (A x ,G y ,C z ) ® (G x ,A y ,F z ). Our measure- 
ments showed that F- and C- modes were at the very limit 
of detection indicating that the Mn moment is essentially 
contained within the a6-plane. Reflections from these 
modes were not included in the final refinements. |25] 
Our analysis yielded a magnetic structure for Mn given 
by the moment in u B : m Mn = Ti (2.83(12), 0.51(4), 0) <X> 
r 3 (0.55(4), 3.79(7), 0). The phase between the two irreps 
is also a variable parameter found to be 0.474(12)7r close 
to -| expected for a cycloid. [25] The values of the Mn- 
spins in this commensurate cycloidal structure along the 
a— and b— axis are given in Fig. |3]for the z=0 and z — g 
layers. 

The magnetic structure indicated by this model is 
dominated by the A modes of these two irreps (A x and 
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FIG. 3: (Color online) Results of refinements of the mag- 
netic structure at 8.5K and H\\b—5T. In panels (a) and (b) 
we show the variation of the M x and M y components of the 
Mn-spins for the ions located at z=0 (black) and z = \ (red) 
along the magnetic propagation vector. Note that the compo- 
nents along a and b are out of phase so as to yield a cycloidal 
ordering shown in Fig. [l] In panels (c) an (d) we show the 
comparison between observed and calculated magnetic struc- 
ture factors (F 2 ) for the analysis of the Mn r = i reflections 
and nuclear plus Tb-magnetic reflections respectively. Since 
the Tb magnetic propagation vector is r =0 the nuclear and 
magnetic reflections are not separated in reciprocal space as 
in the case for the Mn-ions. 



A y ) producing an elliptical cycloid, shown in Fig. [ijb). 
Here the Mn cycloid is contained within the a6-planc 
and Mn-spins rotate around the c— axis, consistent with 
the direction of the ferroelectric polarization along the 
a— axis (P ~ b x c). The anisotropy of the Mn cycloid 
in both low field be and high field ab configurations ap- 
pears to be very similar. In both cases the A y mode pos- 
sesses the higher moment (3.9 and 3.79 hb/Mti respec- 
tively) while the components orthogonal to this mode are 
smaller and of the same magnitude (2.8/xs/Mn). This 
result shows that the flop of the cycloidal plane does not 
effect the anisotropy of the Mn ordering itself. The devi- 
ation of the phase shift between the two irreps away from 
the ideal value produces an angle between spins of 85°. 
Despite this, Si x Sj + i remains parallel to the c— axis and 
should not influence the magnitude of P s significantly. 
Finally for the Mn ordering we find that the G x and G y 
modes are active with amplitudes of ~0.5/ifl/Mn and its 
net effect on the over all cycloidal order that produces a 
polarization along the a— axis is also relatively small. 

We now turn our attention to the Tb-spin ordering. 
We find that in the P\\a phase Tb-spins order com- 
mensurately with the underlying primitive lattice (i.e. 



r-Tb 



0). Analysis of the commensurate Pbnm re- 
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flections clearly showed additional intensity that can be 
modelled by the Tb magnetic order. The best fit to the 
measured data was obtained for a ferromagnetic align- 
ment of Tb-spins along the &-axis and antiferromagnetic 
coupling between nearest-neighbor Tb-spins along the 
[110] direction (Fig. [ijb)). Analysis of this structure 
yields a total Tb moment of 7.24(7) /Ub, with an anti- 
ferromagnetic component of 6.07(9)/is along the a-axis 
and a ferromagnetic component of 3.92(6)/!^ along the 
6-axis. The ferromagnetic ordering along the 6-axis is in- 
deed evident in magnetization measurements under sim- 
ilar conditions [2 r. 

The work we present here unambiguously proves that 
the commensurate P\\a phase in TbMnC>3 coincides with 
an ab Mn spin cycloid for _ff||o=5T. The antisymmet- 
ric DM interaction in this case does yield a ferroelec- 
tric polarization along the a— axis as indeed is observed 
(P s ||a = r x (S, x S»+i) = cxb). We may thus identify 
the inverse DM interaction as the main mechanism for 
the magnetic-field induced flop of ferroelectric polariza- 
tion. 

In a perfect cycloidal magnetic arrangement the ex- 
change mechanism proposed in references [T3J [T3] does 
not yield any ferroelectric polarisation, as the scalar 
product (Si • Sj+i) of neighboring spins is everywhere the 
same. This still holds for the commensurate perfectly 
circular spiral. In the case of an elliptical commensu- 
rate cycloid the exchange mechanism does cause ferro- 
electric polarization as the scalar product varies along 
the modulation. In our case with a modulation of four 
orthorhombic lattice distances the mechanism of [T3J E] 
may thus yield a finite polarisation along the a-direction 
which, however, should still be small due to the only 
minor deviation from a perfect circular cycloid. The ex- 
change mechanism may gain further importance in the 
case of a very anharmonic cycloid. In the extreme an- 
harmonic arrangement, where spins point either along 
the b or a directions yielding the sequence along the b- 
direction shown in fig. [TJc) : left, up, right, down, there 
will be a very effective exchange striction mechanism as 
the scalar product Sj • Sj+i, between nearest neighbor- 
ing Mn ions alternates along the o-axis (in the example of 
fig.[TJc) the alternating angles are 90 and deg. along the 
6-direction) . In the current analysis we have performed, 
the phase difference between rows of cycloids propagating 
along the &-axis (e.g. in Fig. [TJc), rows starting at ions 1 
and 2) is fixed by symmetry as uj — ttt/2 — 45° [23] [27]. 
A deviation from this value will yield an anharmonic spi- 
ral, however such as case can not be fully tested with 
the current data and we can not exclude a small degree 
of anharmonicity. In view of the current debate about 
the electromagnon it appears interesting to add that this 
exchange-striction mechanism will always yield ferroelec- 
tric polarization along the a-direction. 

In conclusion our work demonstrates that the flop in 
the ferroelectric polarization observed in TbMnC>3 arises 



from the flop of the Mn cycloidal plane from be to ab. 
The cycloidal magnetic structure we establish here for 
the high field P\\a phase does not posses a dominant fer- 
roelectric mechanism based on exchange striction. The 
ordering of Tb-spins in this high field phase is that of a 
canted antiferromagnet giving a significant ferromagnetic 
component along the 6-axis. 
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